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Abstract 

Direct CP violation in the h process is a sensitive probe of physics beyond the Standard 

Modeh We report a measurement of the CP asymmetry in B ^ Xs'j, where the hadronic recoil 
system Xg is reconstructed using a pseudo-reconstruction technique. In this approach there is 
negligible contamination from 6 — > ^7 decays, which are expected to have a much larger CP 
asymmetry. We find Acp = 0.002 it 0.050(stat.) it 0.030(syst.) for B Xsj events having recoil 
mass smaller than 2.1 GeV/c^. The analysis is based on a data sample of 140 fb~^ recorded at the 
T(4S') resonance with the Belle detector at the KEKB e+e" storage ring. 

PACS numbers: ll.30.Er, 13.20.He, 14.40.Nd 
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Radiative B decays, which proceed mainly through the b ^ s'j process, have played an 
important role in the search for physics beyond the Standard Model (SM). The inclusive 
branching fraction has been measured by CLEO 0|, ALEPH j^] and Belle jsj, giving results 
consistent with the recent theoretical prediction of (3.57 ±0.30) x 10~^ 0]. The SM predicts 
very small direct CP violation in 6 — 57; the CP asymmetry, 

. , T(b 57) - T(b 57) 

is about +0.5% in the SM, while some new physics, such as supersymmetry, allows the CP 
asymmetry to be above 10% without changing the inclusive branching fraction Thus, 
measurement of Acp may provide information on new physics that cannot be extracted 
from the measurement of the branching fraction. Previously, CLEO measured —0.27 < 
0.965Acp{b 57) + 0.02Acp{b d-f) < 0.10 at 90% confidence level @; however b ^ d-f 
is expected to cancel the CP asymmetry in 6 — > 57 0- The measurement presented here 
has negligible b ^ d'j contamination. 

In this paper, we report on a measurement of the CP asymmetry in i? -^s7, where the 
hadronic recoil system Xg is reconstructed using a pseudo-reconstruction technique. The 
analysis is based on 140 fb~^ of data taken at the T(4S') resonance (on-resonance) and 
15 fb"^ at an energy 60 MeV below the resonance (off-resonance), which was recorded by 
the Belle detector [8] at the KEKB asymmetric e^e~ collider (3.5 GeV on 8 GeV) |9j]. The 
on-resonance data correspond to 152 million BB events. The Belle detector has a three- 
layer silicon vertex detector, a 50-layer central drift chamber (CDC), an array of aerogel 
Cherenkov counters (ACC), time-of-flight scintillation counters (TOP) and an electromag- 
netic calorimeter of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An instrumented iron flux-return for K^/n detection 
is located outside the coil. 

We form the hadronic recoil system, Xg, with a mass up to 2.1 GeV/c^, by combining one 
charged or neutral kaon with one to four pions, where at most one pion can be neutral. We 
also reconstruct Xg via K'^K'^K'^iTr'^) and K^K'^K^^n'^) including the B — > Kcp'j decays 
that were observed recently by Belle |lOj . 

Each of the primary charged tracks is required to have a momentum in the e+e~ center- 
of-mass (CM) frame that is greater than 100 MeV/c, and to have an impact parameter 
within ±5 cm of the interaction point along the positron beam axis and within 0.5 cm in 
the transverse plane. These tracks are identified as pion or kaon candidates according to a 
likelihood ratio based on the light yield in the ACC, the TOP information and the specific 
ionization measurements in the CDC. Por the selection applied on the likelihood ratio, 
we obtain an efficiency (pion misidentification probability) of 90% (7%) for charged kaon 
candidates, and an efficiency (kaon misidentification probability) of 93% (7%) for charged 
pion candidates. Tracks that are identified as an electron or muon are excluded. 

Kg candidates are formed from vr+vr" combinations whose invariant mass is within 
8 MeV/c^ of the nominal K^g mass. The two pions are required to have a common vertex 
that is displaced from the interaction point. The Kg momentum direction is also required 
to be consistent with the Kg fiight direction. Neutral pion candidates are formed from pairs 
of photons that have an invariant mass within 16 MeV/c^ of the nominal 7r° mass and an 
opening angle smaller than 60°. Each photon is required to have an energy greater than 
50 MeV. A mass constrained fit is then performed to obtain the ir^ momentum. 

The B meson candidates are reconstructed from the Xg system and the highest energy 
photon with a CM energy between 1.8 and 3.4 GeV within the acceptance of the barrel ECL 



4 



(33° < 6j < 128°, where 6.y is the polar angle of the photon in the laboratory frame). In 
order to reduce the background from decays of vr'' and r] mesons, we combine the photon 
candidate with all other photons in the event and reject the event if the invariant mass of 
any pair is within 18 MeV/c^ (32 MeV/c^) of the nominal vr*^ (77) mass. 

We use two independent kinematic variables for the B reconstruction: the beam-energy 

constrained mass M^, = ^ {E^^U c^f - {\rxAP*\/ ^? and = E*^^ + E; - E^,,^, 
where -Ebcam is the beam energy, and p*, E*, px^, E*x^ are the momenta and energies 
of the photon and the Xg system, respectively, calculated in the CM frame. In the Mbc 
calculation, the photon momentum is rescaled so that \p*\ = (-Ebeam ~ -^xj/"^ satisfied; 
this improves the Mbc resolution to 2.9 MeV/c^. We require Mbc > 5.24 GeV/c^ and 
-150 MeV < < 80 MeV. We define the signal region to be Mbc > 5.27 GeV/c^. 

The largest source of background originates from continuum e'^e" QQ {q = u, d, s, c) 
production including contributions from initial state radiation {e~^e~ QQ'j)- In order to 
suppress this background, we require the presence of a high energy lepton from the opposite 
B. The lepton can be either an electron with a CM momentum greater than 0.8 GeV/c or 
a muon with a laboratory momentum greater than 1.0 GeV/c. In both cases, the opening 
angle 6^^ between the high energy photon from the signal B meson and the lepton, calculated 
in the CM frame, must satisfy | cos 6*1^1 < 0.8. In addition to this lepton requirement, we 
use the likelihood ratio described in Ref. ^3], which utilizes the information from a Fisher 
discriminant l2| formed from six modified Fox- Wolfram moments 3] and the cosine of the 
angle between the B meson flight direction and the beam axis. The lepton requirement 
(likelihood ratio selection) retains 15% (92%) of the signal, rejecting 98.5% (55%) of the qq 
background. 

When multiple B candidates are found in the same event, we take the candidate that 
gives the highest confidence level when we fit the Xg decay vertex, constrained to the profile 
of the measured interaction region. The confidence level for the Kg^^'j mode is set to 
zero, because we do not determine the vertex. If multiple B candidates appear in an event 
according to the inclusion or omission of different 7r° mesons in the X^ recoil system, we 
take the candidate that has the minimum |Ai?|. 

We fit the M^c distribution to extract the signal yield. The Mbc distribution of the qq 
background is modeled by an ARGUS function |lj| whose shape is determined from the off- 
resonance data. (The lepton requirement is not applied to the off-resonance data in order 
to compensate for the limited amount of data in that sample.) Background from B decay is 
divided into two components. B decays through b ^ c transitions (except color-suppressed 
B decays such as 5° — >■ D^ir^), which we call BB background in this paper, have a non- 
peaking Mbc distribution which is modeled by another ARGUS function. Rare B decays, i.e. 
B decays through b ^ s and b u transitions (charmless B decay) and color-suppressed B 
decays, are not negligible and have peaks at the B mass in the Mbc distribution. The sum 
of these distributions is modeled by a Gaussian plus an ARGUS function. The shape of the 
distributions is determined by a corresponding Monte Carlo (MC) sample. 

The Mbc distribution of the signal component is also modeled by the sum of a Gaussian 
and an ARGUS function. All parameters are determined from the B Xgj signal MC 
simulation, except that the mean of the Gaussian is extracted from a fit to i? ^ Dtt decays, 
described below. Our nominal signal MC contains a i? — *• K*{892)'j component and an 
inclusive b s'-f component with Mx, > 1.15 GeV/c^. The ratio of the two components 
is based on the branching fraction for B f^* (892)7 measured by Belle and the PDG 
value of B{B Xg'j) [16]. The Xg system of 6 — > 37 is modeled as an equal mixture of sd 
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and su quark pairs, and is hadronized using JETSET [17|. The Mx^ spectrum is fitted to 
the model by Kagan and Neubert 0] with the b quark mass parameter nif, = 4.75 GeV/c^. 

Before discussing the extraction of Acp, we elaborate on the signal modeling. We create 
alternative signal MC samples in which the i^* (892)7 fraction is varied by ±la or the b 
quark mass parameter is varied by ±0.15 GeV/c^. We account for the uncertainty in the 
hadronization process by preparing an alternate MC sample wherein B candidates with 
> 1.15 GeV/c^ are selected or discarded at random to match the fractions of such 
candidates observed in the data having between one and four pions. We prepare other such 
samples wherein the proportion of selected B candidates without a neutral kaon or pion 
matched the value seen in data; the correction is modest for all modes except Kn'-f (which 
contributes 32% in the MC but 12% after correction). The systematic error estimate based 
on these samples is described below. 

In the pseudo-reconstruction analysis of B —>■ Xsj, the flavor of the B meson is "self- 
tagged" when the net charge of the Xg system is non-zero or the Xg system contains an 
odd number of charged kaons. Otherwise the flavor is "ambiguous." Because we have three 
possibilities for the flavor tag, there are three ways to tag the flavor incorrectly: Wi is the 
probability of classifying a self-tagged event as a self-tagged event of the opposite flavor; W2 is 
the probability to classify an ambiguous event as self-tagged; ws is the probability to classify 
a taggable event as ambiguous. For example, W2 refers to the case when a B^ ^ Kgir^TT^'j 
event is tagged as a Kgn'^ir^'y event, and refers to the case when a B~^ Kgn'^ir'^'y event 
is tagged as a Kg7T~^7r~'y. The effect of the flavor dependence of the wrong tag fractions is 
negligibly small, and is included in the systematic error from the B Dn study described 
later. 

The formula to calculate Acp is then Aqp = -DA^p with the dilution factor D = (1 — 
W2 - W3)/[{1 - W2)(l - 2wi - ws)] and the raw asymmetry A^^p = {N_ - N+)/[N^ + N+ - 
{W2/ (1 — W2))Nq\, where iV_ (A^+) is the number of events tagged as originating from a b (b) 
quark, and A^o is the number of events classified as ambiguous. 

The three wrong tag fractions are estimated using signal MC, and thus are model de- 
pendent. We calculate the wrong tag fractions for the collection of signal MC samples 
described earlier, and extract the systematic errors in the wrong tag fractions from the 
changes in Wi as the model is changed. We find Wi = 0.0206 ± 0.0027, W2 = 0.248 ± 0.020 
and ws = 0.0067 ± 0.0013, resulting in a dilution factor of D = 1.041 ± 0.006. 

Figure HI shows the Mbc distributions for events that are classified as 6-tagged, 6-tagged 
and ambiguous, respectively. The distributions are fitted with signal, qq, BB and rare B 
components. We assume the shapes of the three components are common for b- and 6-tagged 
classes, distinct from those for the ambiguous class. In the fit, the numbers of events from 
the BB and rare B background in the signal region are fixed to be 39.3 (9.2) and 35.4 
(3.3) for the common b- and 6-tagged classes (ambiguous class) using the MC prediction, 
while the normalization of the qq component is allow to float. Signal yields are obtained 
by integrating the signal Gaussian and ARGUS functions in the signal region. We find the 
signal yield in each class to be = 393.2 ± 25.9, = 392.0 ± 25.9 and Nq = 52.8 ± 9.6, 
resulting in Acp = 0.002 ± 0.050. 

The systematic errors are studied using a. B —>■ Dn control sample. D mesons are 
reconstructed in the same final states as Xs except for the final states to which D cannot 
decay. Only a mild requirement on the invariant mass of the reconstructed D mesons of less 
than 1.9 GeV/c^ is applied, in order to allow cross-feeds between different D decay channels. 
In the calculation of Mbc, the primary pion momentum is rescaled in the same way as the 
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FIG. 1: Mbc distributions for (a) events tagged as b, (b) events tagged as b and (c) events classified 
as ambiguous. Fit results are overlaid. 



photon in the B Xg^ reconstruction. 

In order to check the vahdity of the signal Mbc shape obtained from MC, we compare the 
Mbc distribution for B — > Dtt data and MC. The data Mbc distribution is fitted with the 
sum of a Gaussian and an ARGUS function plus qq and BB distributions. The qq and BB 
distributions are obtained from off-resonance data and MC simulation, respectively, where 
each normalization is scaled according to the luminosity. Prom this procedure, we determine 
the mean of the Gaussian used to model the signal Mbc distribution. We find that the ratio 
of the Gaussian and ARGUS function agrees between data and MC within the statistical 
error of 23%. 

The systematic error on Acp due to the fitting procedure is estimated by varying the 
value of each fixed parameter by ilo" and extracting new signal yields and Acp for each case. 
We assign an additional 23% error obtained from B Dn to the ratio of the Gaussian and 
ARGUS function. The ratio is also varied separately for and to take into account 
the possible difference of sub-decay modes between b and b. We also use the signal shape 
parameters obtained from each of the signal MC samples described earlier and extract Acp 
for each case. We set the normalization of either the BB or rare B backgrounds to zero 
and to twice its nominal value to account for its uncertainty. A 50% error is assigned to the 
Gaussian width and to the ratio of the Gaussian and ARGUS function for the rare B decays, 
to compensate for our limited knowledge of their branching fractions. The changes of Acp 
for each procedure are added in quadrature, and are regarded as the systematic error. We 
obtain a systematic error of 0.014 due to the fitting procedure. 

The B Dn sample is also used to estimate the possible detector bias in Acp- We 
estimate Acp for B — > Dir using the same method as for B — > Xg'y. The flavor is determined 
from that of the D candidate. Prom the MC study, we find the B — > Dn wrong tag fractions 
to be wi = 0.0314 ± 0.0039, W2 = 0.285 ± 0.024, and = 0.0225 ± 0.0033, resulting in a 
dilution factor of D = 1.059 ± 0.013. The signal yields are calculated by fitting the data 
Mbc distribution with signal, qq, and BB components with fixed BB normalization, and 
are found to be = 2125 ± 60, N+ = 2105 ± 60 and A^o = 580 ± 30, resulting in a CP 
"asymmetry" of Acp — 0.006 ±0.022. We therefore assign 0.022 as the systematic error due 
to detector bias. 

In addition, we estimate the uncertainty due to possible asymmetries in charmless B 
decays. We divide these decays into two groups: those for which Acp has been measured 
{B — > K*rj, Kn^ and Kp), corresponding to about half of the selected charmless B events; 
and all other decay modes. We assign an Acp of ±30% to the first group (corresponding to 
around la in each case), and ±100% to the second, assuming the signs oiAcp are correlated. 
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FIG. 2: Aqp for each Mx^ bin. Mx,, is the reconstructed value of the mass of the Xs system. The 
data points do not include any correction for the effect of smearing and mis-reconstruction. 



The resulting contribution to the measured Aqp, 0.014, is taken as a systematic error. 

The systematic error on Acp due to the uncertainty of the wrong tag fractions is found 
to be small. The contribution of the b ^ dj process is negligible, because we require the 
existence of kaons in the final state. This is confirmed by examining MC samples of -B — > p7, 
LU'j and inclusive b d'j process we expect that the contribution of b ^ d'j to Acp is 
less than 0.001. The total systematic error on Acp is then calculated to be 0.030 by adding 
in quadrature the errors mentioned earlier. 

In order to examine the Mx^ dependence of Acp, we divide the data sample into six bins 
of measured Mx^ and perform the Mbc fit for each bin. Figure |21 shows the Aqp distribution 
as a function of Mx^- Here, the systematic error due to the detector bias and possible 
CP asymmetry in charmless B decays is not included in the error. We find that Acp is 
consistent with zero for all Mx^ values in the distribution. From our MC study, about 2% 
of the events reconstructed with Mx^ < 2.1 GeV/c^ are expected to have a true Xg mass 
greater than 2.1 GeV/c^. 



In conclusion, the CP asymmetry of B 
measured to be 



for events with Mx, < 2.1 GeV/c^ 



IS 



Acp{B XsT, Mx, < 2.1 GeV/c^) = 0.002 ± 0.050(stat.) ± 0.030(syst.), 

consistent with no asymmetry. This corresponds to —0.093 < Acp < 0.096 at the 90% 
confidence level, where we add the statistical and systematic errors in quadrature and assume 
Gaussian errors. The result can restrict the parameter space of new physics models that 
allow sizable CP asymmetry in 6 — >■ 57 j2ol |. 
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